Evolution of the Fermi surface in phase fluctuating d-wave superconductors 
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One of the most puzzling aspects of the high T c superconductors is the appearance of Fermi arcs 
in the normal state of the underdoped cuprate materials. These are loci of low energy excitations 
covering part of the fermi surface, that suddenly appear above T c instead of the nodal quasiparticles. 
Based on a semiclassical theory, we argue that partial Fermi surfaces arise naturally in a d-wave 
superconductor that is destroyed by thermal phase fluctuations. Specifically, we show that the 
electron spectral function develops a square root singularity at low frequencies for wave-vectors 
positioned on the bare Fermi surface. We predict a temperature dependence of the arc length that 
can partially account for results of recent angle resolved photo emission (ARPES) experiments. 
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The transition from a superconductor to a normal 
state at T c can take place in two qualitatively differ- 
ent ways. Conventional metallic superconductors take 
the BCS-Eliashberg route by which the amplitude of the 
complex order parameter |A|e Z¥3 vanishes at T c together 
with the quasi-particle gap. Thus the transition is accom- 
panied by formation of a normal Fermi surface. Phase 
fluctuations play a negligible role in this transition and 
T c is determined by the zero temperature gap Ao. On 
the other hand, in superconductors with small superfiuid 
density n s , the phase-stiffness defines a temperature scale 
Tg ~ n s /m* that may be much lower than A. The tran- 
sition then occurs at T c ~ Tg due to disordering of the 
phase by thermal fluctuations. This is thought to be 
the case in the underdoped cuprate superconductor*^, 
where T c is seen to increase with doping concomitantly 
with the superfiuid density^ and is anti correlated with 
the quasi-particle gap. Is such a transition, driven by 
phase fluctuations, accompanied by formation of a Fermi 
surface as is the BCS case? 

Naively, the answer is negative. At temperatures not 
too high above T c (T < Tmf ~ Ao) pairing survives in 
the normal state. One might expect that in this regime 
the single particle spectral function remains peaked at 
finite frequency of order A , the energy needed to dis- 
sociate a pair. This reasoning is certainly correct for 
an s-wave superconductor with clear separation of scales 
between the phase stiffness and the larger quasi-particle 
gap. Does it still hold in a d-wave superconductor, where 
such a clear separation of scales cannot exist because of 
the vanishing gap at the nodal points? 

In this paper we investigate the evolution of the Fermi 
surface in the normal state of a d-wave superconductor by 
studying the coupling between nodal quasi-particles and 
strong phase fluctuations. Our main result is that for any 
value of the anti-nodal gap a sharp Fermi surface forms 
in the normal state immediately above T c . Nevertheless 
this is not a normal Fermi liquid. Because of coupling 
to thermal phase fluctuations, the quasi-particle peak on 
the Fermi surface is replaced by a square root singularity 
in the spectral function A(k = k/,w) = a(k)/ v / w. This 
singularity is present everywhere on the Fermi surface, 
but at wave-vectors far from the gap nodes, that is 5k > 



\J (V<y5 2 ) = qo(T), it is exponentially suppressed by the 
pre-factor a(k). Therefore a peak at zero frequency is 
observable in practice only along arcs of length ~ go on 
the Fermi surface, centered at the nodal points. The 
dominant spectral feature further from the node, is a 
broad peak at frequency uj = A^. 

Broadening of the quasi-particle peak at u> = Ak by 
thermal phase fluctuations in the normal state was noted 
in previous analyses^. However the appearance of a 
singularity at zero frequency signaling the formation of 
a Fermi surface, was overlooked. We argue that this 
physics underlies the formation of the Fermi arcs seen 
in angle resolved photoemission spectroscopy (ARPES) 
of underdoped cuprates^. 

Before proceeding we note that a number of recent 
works proposed alternative mechanisms for the forma- 
tion of Fermi arcs and their evolution with temperature. 
Varma and Zhu£ have argued that current fluctuations in 
a state with broken time reversal symmetry can explain 
the observed phenomena. Paramekanti and Zhao^ postu- 
lated that the normal state of the underdoped cuprates is 
dominated by a quantum critical point with deconfmcd 
spinons at the critical doping level for onset of supercon- 
ductivity. In their theory, the broadening of the spectral 
peaks accompanying the deconfinement leads to forma- 
tion of the arcs. Kim et al^ have discussed the formation 
of Fermi arcs due to critical fluctuations of an underly- 
ing nematic order parameter in the presence of quenched 
disorder. 

For concreteness and to make contact with experiments 
in the cuprates, we consider a layered structure relevant 
to these materials. The main ingredient in our analy- 
sis, leading to emergence of a Fermi surface, is the cou- 
pling of the dirac quasi-particles to the local currents 
j s = ?i s V(/j induced by the phase fluctuations. These cur- 
rents are conveniently split into two contributions with 
markedly different effect on the quasiparticles: A trans- 
verse part originating from vortex configurations and an 
irrotational, or longitudinal part. We find that the in- 
teresting effects on the single particle spectral function 
result from the transverse fluctuations. To simplify the 
subsequent analysis we treat the vortices as purely ther- 
mal fluctuations and neglect their dynamics altogether. 
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FIG. 1: In the presence of static current, the dirac cone, seen 
in panel (a) on the bare Fermi surface, is shifted in energy 
and boosted by the wave- vector q associated with the super- 
current. If the current is parallel to v/ as in (b,c) the net 
effect is that the cone slides up or down leaving one point 
of zero energy excitations invariant. If qlv/ as in (d) the 
nodal point is transported along the bare Fermi surface. As 
shown in panel (a) there is a continuous range of currents q 
that lead to a zero energy peak on each point of the bare 
Fermi surface. This pile up of spectral weight results in a 
singular peak at zero energy on the Fermi surface. 

This is a reasonable assumption in the normal state where 
vortex motion is expected to be diffusive and slow on elec- 
tronic time scales. The longitudinal fluctuations, by con- 
trast, can be shown to be highly quantum 11 . However, 
we will show that their coupling to the quasi-particlcs is 
irrelevant, even if the small gap in their dispersion, due 
to the c-axis Josephson plasmon, is neglected. 

Under these assumptions, the electron Green's func- 
tion relevant to ARPES is given by: 

G<<T,t) = J Pj ± (r)P[u(r)] (^(0, 0)^(r, t)) {Ji(r)} . 

(1) 

The external integral implements a thermal average over 
the transverse currents (or vortex configurations). The 
thermal weight of a current configuration, P[jj_(r)], is 
determined by a classical x — y modeli 2 -. The internal 
average is the Green's function in a given static current 
configuration. This includes integration over the dynam- 
ical longitudinal phase fluctuations. 

In the spirit of refs. [HQ we will approximate (JTJ) by 
an average over a distribution of uniform currents. This 
amounts to a semiclassical approximation valid if spa- 
tial variations of the current over a Fermi wavelength are 
small (i.e. kf£ 1). Interference effects due to coher- 
ent scattering of quasi-particles by vortices are neglected 
in this scheme, which is justified at the temperatures of 
interest in the normal state. 

With these approximations, the spectral function on 
the two dimensional plane is given by: 

AQt,u) = f dV(q)A q (k,^). (2) 



Here ^4 q (k, to) is the spectral function in the presence of a 
uniform supcrcurrent parameterized by a wave- vector q. 
Such current corresponds to a pairing field that depends 
on the center of mass coordinate of the pair as A k e* 2qR . 

We may now use the standard BCS formalism in order 
to calculate A q (k, tu) for a given supercurrent q. This is 
easily seen to be 

A q (k,w) = |u k _ q | 2 (5(cj-£; k _ q (q))+|i; k „ q | 2 ( 5(tJ+P q _ k (q)) 

(3) 

and the quasi-particle dispersion to leading order in q is 
given by 



£ k (q)«v / (k)-q+ v /£2 +A 2. (4) 

for k near the bare Fermi surface, defined by £ k/ = e k/ — 
fi = 0. To this order in q, the coherence factors u k 
and v k are unchanged by the current, that is |it k | = 
(1 + UV& + ^L)/^ and K| 2 = 1 - | Mk | 2 . 

The essential physics leading to emergence of a Fermi 
surface in the normal state can be understood qualita- 
tively from Eqs. j3]) and ((4]). Fig. [T] shows the dispersion 
of the quasi-particle peaks in A q (k, u) for different direc- 
tions of the supercurrent. The momentum of the peak 
is shifted by q because of the boost associated with the 
supercurrent, but its frequency is also "Doppler shifted" 
by V/ • q. For a current parallel (or anti-parralel) to V/, 
this amounts to the dirac cone sliding up (or down) on 
the slope of the original dispersion as seen in Fig. []Jb,c). 
The important thing to note is that in both cases there 
is a zero frequency peak, or "Fermi point", whose posi- 
tion is independent of the current. On the other hand, 
for current perpendicular to vj there is no doppler shift, 
and the node is simply transported a distance q along 
the Fermi surface (Fig. [TJi). In general an arbitrary 
current q will always lead to a zero frequency peak lo- 
cated precisely on the bare Fermi surface. The parallel 
component of q leaves the zero frequency peak in place 
on the nodal point, while the perpendicular component 
transports it on the fermi surface. In fact, as shown in 
Fig. QJi, for every point on the Fermi surface, there is a 
continuous range of wave-vectors q that lead to a zero 
frequency peak at that point. This pile up is the ori- 
gin of the zero frequency singularity that develops on the 
Fermi surface. However, the probability for a current 
fluctuation with wavevector larger than a typical qo is 
strongly suppressed in the distribution -P(q). Therefore 
on points of the fermi surface at a distance larger than 
qo from the node the singularity is multiplied by a very 
small pre-factor making it essentially unobservable. 

We now turn to a systematic derivation of the spec- 
tral function using Eqs. ©,© and We concen- 
trate on wavevectors k sufficiently close to the nodes so 
that we may approximate the dispersion by a Dirac cone. 
Namely: A k = va ■ (k— k„) and £ k = v/ ■ (k— k n ). Let us 
also change variables from q to rj = Vf ■ q and C = va • q. 
Plugging these definitions to ^ we have 
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where for notational simplicity we omitted the subscript 
k from Ak and £k- The distribution P(rj,() is obtained 
by a simple change of variables from -P(q), which is taken 

to be a gaussian of the form P(q) = A/" exp(— q ^ T ^ ) [|| . 
We will later discuss the temperature dependence of ga- 
in Ref. Q, the integral (fT5|) was approximated by 
expanding to leading order in 77/ Ak and C/Ak- For k on 
the Fermi surface far from a gap node this approximation 
does give the the dominant spectral feature, which in this 
region is simply a broadened peak at frequency uj ~ Ak. 
However, the expansion cannot be generally valid for a 
superconductor with gap nodes, on which these ratios 
diverge. In particular, it cannot capture the pile up of 
low energy spectral weight on the Fermi surface due to 
movement of the Dirac cone in energy and momentum, 
as illustrated in Fig. [TJ 

We proceed to evaluate (fT5)) without resorting to an 
expansion in rj/A^ and C/Ak- In what follows we out- 
line the main steps in the derivation and relegate the 
details to the online supplementary information. First 
we transform to the polar coordinates (r, 9) defined by 
r = v / (^-^) 2 + (A-C) 2 and cosfl = (f - r/)/r. This 
allows to resolve the delta functions in (|15|) by integrat- 
ing over r. Next to facilitate the angular integration we 
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FIG. 2: Calculated ARPES energy distribution curves (EDC) 
for different wave- vectors k on the Fermi surface at increasing 
distance from the node. The curves are obtained by numerical 
integration of Eq. (|20|l , multiplied by the Fermi function and 
convolved with a narrow gaussian function to mimic the finite 
experimental resolution. The dashed line follows the disper- 
sion of the peak, demonstrating the formation of a Fermi arc. 
Evidently the peak departs from ui = only at k — k n w 1.5(70 • 



make another change of variables via a = s cos 9. Now, 
for wavevectors k on the Fermi surface (i.e. £ = 0) the 
integral reduces to a rather simple form: 



M 
2 



1 (1 -a) 1 / 2 
da [ n ; a {, /9 F(\u;\ : a) (6) 

1 



■(1 + a) 5 / 2 " 



where 



F(\u>\,a) = Afexp 



M 2 a 2 



A 



v 2 q 2 {l + af 



9 9 

V A% 



(7) 



and AT = l/{-KVfV&q 2 ). For \ui\ ^ 0, the divergence at 
a = — 1 in the integrand of ([20]) is cut off by the first term 
in the exponent in ([7]) , which effectively introduces a new 
lower bound to the integral. The singular contribution 
to the spectral function is now easily evaluated as: 



,(k,o;)=A ( UAg o)- 1 e-^/^y» 



(8) 



which is the advertised zero frequency singularity. Aq is a 
constant of order unity. Note that the singular contribu- 
tion comes from a = — 1, which in the original variables 
corresponds to the range of currents q defined by the 
condition q • va = Ak- This is precisely the condition 
for wave-vectors q illustrated in Fig. that pile up low 
frequency spectral weight on the point k on the Fermi 
surface. If k is far from the node compared to a typi- 
cal current fluctuation qo, then A/vAqo is large and the 
singularity is strongly suppressed by the exponential pre- 
factor in pTj) . Then the dominant spectral feature is a 
broad peak at frequency uj ~ A obtained from taking the 
saddle point of l[2Tj]). 

The ARPES spectrum is related to the spectral func- 
tion via /(k, uj) = nj?(a;)^4(k, ui). In Fig. [2] we plot this 
function for k points on the Fermi surface at varying dis- 
tance from the node. The curves are obtained by numeri- 
cal integration of (|20[) and the convolution with a narrow 
gaussian function which mimics the experimental resolu- 
tion. The overall picture is appealingly similar to what 
is seen in ARPES experiments in the normal state of un- 
der doped cuprates^. It is seen that the sharp peaks are 
pinned to zero frequency along a line on the Fermi surface 
up to a critical distance from the node, which is of the 
order of qo. At larger distance from the node the peak 
departs from zero frequency, quickly approaches u = Ak 
and broadens considerably 

A similar analysis as outlined above can be carried out 
for points k on the line of nodes Ak = away from the 
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Fermi surface. In this case the singular contribution to 
the spectral function is given by 



,(k,w) = A (v A q y 



"/go 



(9) 



which appears as a sharp peak dispersing as £k- This too 
is in qualitative agreement with ARPES. 

We turn to discuss the nature of the phase field in 
the normal state and the contribution of different types 
of fluctuations to the electron spectral function. The 
phase field of a fluctuating 2d superconductor can be 
decomposed into a longitudinal (non singular) part and a 
transverse, or vortex contribution. These components are 
governed by independent actions in the long wavelength 
limit. 

Longitudinal phase fluctuations - These fluctuations 
may be treated at the gaussian level. Physically they 
correspond to Josephson plasmons in the d-wave super- 
conductor, which were studied in detail by Paramckanti 
et al [llT | . The structure of their spectrum in a layered 
system is as follows. At low momenta there is a wide 
regime of linear dispersion, with a small plasmon gap 
due to the c-axis Josephson coupling. Experimen- 
l - 10K in Bi2212i^i^, which is lower 
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than typical temperatures of interest in the normal state 
of lightly undcrdopcd cuprates. we shall therefore ne- 
glect this gap. At the high momentum side, the linear 
plasmon dispersion terminates at a characteristic energy 
scale Up at q ~ In the cuprates u p is a high energy 
scale of order a few eV. Therefore, quantum dynamics 
of the plasmons must be taken into account fully at the 
temperatures of interest. 

How do these fluctuations affect the fermion spectral 
function? The most relevant coupling between the Dirac 
Fermions and the longitudinal phase fluctuations is the 
current-current coupling^ of the form Vf ■ \7tp , ip cr i[> C r. The 
scaling dimension of this coupling is easily seen to be 
— 1/2 making it irrelevant at low energies. Further scaling 
arguments, and direct calculation, show that the life-time 
of quasi-particles due to this coupling diverges at low 
energies as 1/r ~ (e — ep) 2 . At finite temperature, this 
will be cut off at T 2 , which gives thermal broadening of 
the quasi-particles. However a Fermi arc does not form 
due to coupling to the longitudinal phase fluctuations. 

Transverse phase fluctuations- The fluctuations con- 
sidered in the calculation of the electronic Green's func- 
tion (P) are the transverse (vortex) contribution. Vortices 
are macroscopic objects whose motion is expected to be 
ovcrdamped. This picture is supported by the fact that 
the measured Nernst signal and diamagnetism in under- 
doped cuprates^ are consistent with the predictions from 
a classical x — y model with overdamped dynamics^. In 
writing ([1]) as a static average, we assumed this dynamics 
to be sufficiently slow on quasi-particle timescales. 

The arc length at T > T c is of order q (T) = {Vlp\_) 1/2 . 
That is, the width of the of the transverse current distri- 
bution, which we shall calculate within the classical two 
dimensional x — y model. We note that in the equivalent 



coulomb gas model, qo is related to the vortex density 
correlation: q 2 = (47r 2 /f2) (n(q)n(— q))/q 2 - 

This quantity was previously estimated^ in the 
high temperature limit using the Debye-Hiickel 
approximation^. Here we present a different cal- 
culation of qo (T), using a variational approach. Our 
result coincides with the Debye-Hiickel approximation 
at high temperature, but can also be used at lower 
temperatures, closer to Tbkt- First we carry out the 
usual mapping of the coulomb gas onto the Sine-Gordon 
model: 
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(10) 

Here K = J/T is the bare phase stiffness of the origi- 
nal x — y model, y = exp(— E c /T) where E c is the vortex 
core energy, and £ c the core size. \x v (x) is a local chemical 
potential that couples to the vortices, and allows the cal- 
culation of vortex density correlations. We shall now ap- 
ply the self consistent harmonic approximation (SCHA) , 
which is known to give good results for the correlations 
in the normal (i.e. gapped) phased. In this approach, 
the cosine term in (|10p is replaced by a quadratic mass 
term \m 2 B 2 . The mass is determined by minimizing the 
variational free energy F$ + T(S — 5*0)0, where Fq is the 
free energy of the quadratic variational action 5*0 and 
()o denotes a thermal average with respect to Sq. Given 



the solution m (T) = j- (SyK) "V 1 ~^~) , it is now straight 
forward to compute qo within the variational action. 



q 2 o(T) 



/Mx|) = 



/oo 
drr In r sinh [2irK f (m£ c r) 

d 2 q e"^ 
| q |< * 2tt q 2 +to 2 ' 



(11) 



Eq. ([36]) can be shown to be of the Debye-Hiickel form 
in the high temperature limit. 

The parameters J, E c , and £ c that control the tempera- 
ture dependence of go, are not easily connected to observ- 
able properties of the cuprates. The parameter J is the 
bare superfluid stiffness of the x — y model. Nevertheless, 
in a d-wave superconductor it is expected to be temper- 
ature dependant J(T) = J(0) — AT due to physics that 
lies out side of the pure x — y model, namely depletion of 
the condensate by quiasiparticlcs at the gap nodes^ -. In 
the normal phase, the temperature dependence may well 
be more complicated because of the emergent finite den- 
sity of states at zero energy. However, to leading order in 
the current fluctuations we may assume that the linear 
decrease of J persists in the normal state. We take the 
parameter A from experiments, that measure the lead- 
ing temperature dependence of the superfluid density^ 2 -. 
Note the distinction between J(T) which is the bare stiff- 
ness of the x — y model and the macroscopic stiffness, 
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FIG. 3: The Fermi arc length is related to the width qo of the 
current distribution. This is plotted as a function of temper- 
ature using eq. (f36|) . for several values of E c /Tbkt (see text 
for details). 

which vanishes in the normal state due to proliferation 
of vortices. 

The results of the calculation of the arc length ([36)) 
are plotted in Fig. [3] for a range of values of E c /Tbkt- 
We assumed that £ c = 5a. Near Tbkt, go (T) is expo- 
nentially suppressed. At higher temperatures, there is a 
region where go increases rapidly due to the proliferation 
of free vortices, followed by a roughly linear increase re- 
gion. Since we have assumed that J(T) = J(0) — AT, 
the bare stiffness vanishes at some temperature (which is 
about 3Tbkt for our chosen value of A = 0.3). At this 
temperature go (T) diverges logarithmically, and the con- 
tinuum description of the x — y model breaks down. The 
temperature dependence of the bare stiffness J above T c 
in the cuprates is not clear, even though there is some 
evidence that it continues to decay linearly over a con- 
siderable ranged. 

In our treatment, we have neglected the inter-planc 
Josephson coupling. This coupling becomes relevant at 
some T c > Tbkt, where a phase transition to a three 
dimensional ordered state occurs. At this temperature, 
our treatment of a two dimensional x — y model is no 
longer valid, vortex formation is strongly suppressed, and 
we expect go (T) to drop abruptly to zero. Indeed, the 
observed arc length seems to drop to zero at T c Q . 

Discussion and conclusions - Before concluding let us 
remark that, we have so far not payed special attention 
to effects of proximity to the Mott insulating state at half 
filling. The simplest (though not necessarily the correct) 
way to include the effect of the strong local repulsion, 



is via the slave boson mean field theory. Within this 
approach, the current carried by a quasiparticlc is renor- 
malized b y a factor a, which is proportional to the hole 
doping x [20| | . Hence the doppler shift of the quasiparti- 
cle dispersion ^ is renormalized by a similar factor, so 
that v/ • q — > qv f ■ q. This would have a dramatic effect 
on our earlier considerations. In particular, application 
of a static current with q parallel to v/ would no longer 
leave invariant the point of zero energy excitations. Ac- 
cordingly, the average over the current distribution ([3]), 
would not result in a singular peak on the Fermi sur- 
face. Specifically, we find the leading edge of the EDC 
on the Fermi surface behaves as \u>\ in this case. It is 
interesting to note that the same renormalization factor 
a, which leads to the present discrepancy of slave boson 
theory with experiment, is also responsible for its famous 
failurrj^ to explain the doping independent slope of the 
superfluid density with temperature 2 ^. 

To summarize, we have shown that destruction of a 
d-wave superconductor by proliferation of thermal vor- 
tices must be accompanied by formation of at least a 
partial fermi surface. We argued that this phenomenon 
is fundamental to the appearance of Fermi arcs in the 
normal state of underdoped cuprates^, and may partially 
explain their evolution with temperature^. The main re- 
sults rely on two central assumptions. The first is that 
vortices are purely thermal and their quantum dynam- 
ics may be neglected. The second is the calculation the 
fermion spectral function in the presence of static vor- 
tices within a semiclassical approximation along the lines 
of [3] . The semiclassical approximation was tested in nu- 
merical simulations 2 ^, which verified that it is controlled 
by the small parameter v^/vf. It therefore seems safe to 
apply the semiclassical approach to the high T c cupratc 
superconductors, except in the very underdoped regime. 

In extremely underdoped materials quantum dynamics 
of the vortices is also expected to be increasingly im- 
portant, in violation of our first assumption. Indeed, 
the onset of superconductivity possibly corresponds to 
a quantum phase transition, driven by fluctuating quan- 
tum vortice o 24 ' 25 . An important open question is how 
the normal state we describe, that arises in a thermal 
vortex liquid, evolves to the highly quantum regime near 
the critical doping for the onset of superconductivity. 
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Supplementary information 



I. DERIVATION OF THE ANGLE RESOLVED PHOTOEMISSION (ARPES) WEIGHT 

A. Static current 



A static current is parameterized by a wave- vector q that can be ascribed to a twist in the order parameter, or 
to the hopping matrix elements of the Hamiltonian as an external vector potential. The two descriptions are related 
by a simple gauge transformation. The ARPES spectrum in absence of perpendicular magnetic field is given by the 
gauge invariant Green's function 



I(k, u) = i J drdtG< a (r, t) exp (+ik ■ r + iivt - i j A • dl^ 



where G< a (r, t) = {c\ o (0)c a ^(t)) . We induce a current using the gauge A = q, for which 



/ q (k,«)= / ^e^(4 k _ q (0)c CT)k - q (t)) 



(12) 



(13) 



Note that this is simply the equilibrium expression for G < , boosted by momentum q. We express (1131) in terms of 
the Bogoliubov operators a k that diagonalizc the Hamiltonian: 

7 q (k, u) = |u k _ q | 2 (aj c _ q)T a k - q ,-r)£(w - #k- q (q)) + K- q | 2 (a- k+q ,j.al k+q!i )<5(w + £ q _ k (q)) 
= n F (cj) (|u k _ q | 2 <5(>> - £k_ q (q)) + K- q | 2 <5(^ + £ q _ k (q))) 



n F (w)A q (k, -u>). 



(14) 



Here -E k (q) = v/(k) • q + \/£ k + A k is the quasi-particle spectrum to leading order in the current. Note that the 
coherence factors u k and i; k are unchanged by the current to this order in q. That is, |u k | 2 = (1 + £ k /.E k (0))/2 and 



(1 - &/£k(0))/2. 



B. Fluctuating current 



As we discuss in the paper, the spectral function in the presence of the phase fluctuations is obtained by averaging 
l4|) over a distribution of static currents. This leads to the following integral: 



A(k )W ) = i J d V d(p( v x) £ ^i 



V(£ - V) 2 + (A - C) 5 



6(u-ri- s V (^) 2 + (A-C) 2 ), 



(15) 
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where we have denoted: A = va • (k — k„), £ = v/ ■ (k — k n ), rj = v/ • q, and £ = va ■ q- The current distribution is 
assumed to be gaussian, such that 



p(v,0 



>r 



9 ^Xp I 

TTv f v A q£ V vfqfi v A q$ 



(16) 



To compute ([15]) we change to the polar coordinates rcos6> = £ — r), rsinf? = A — Then (|15p takes the form 

A(k,u) = - I d6drrP(. A {r,6) V (1 - s cos 0)5 (w - £ + sr(l + scos6)) (17) 

2 J B=±l 



where P^.a(t, 9) = P(£ — rcosO, A — r sm0). Integrating over r to resolve the delta functions and changing variables 
to a — s cos 6 we have 



1 f\(l-a)V* 



where 



1 



2 exp < 



2 \ 



2 2 



2 2 

«a?o 



(19) 



For wave-vectors k on the Fermi surface we plug £ = into (fT8|) and (fT9|) . In this case fo,— A(k>> <*) = -Po,a( — a) 
and therefore (|T5]) is reduced to 



A(k,o,) = M £d a |L^F 0lA (M l a) 



(20) 



The divergence at a = —1 is cutoff by the first term in the exponent of the distribution -Fo.a(M,oO which leads to 
the singular contribution at zero frequency: 



^ S (k,U,) = ( V A S 0)-V^o^ 



on the other hand, on the nodal line (Ak = 0) l|18p simplifies to 

",J1^) 1/2 
a) 5 / 2 



1 f (1 

AQl,w) = -\w-Z\ J ^da [ — 



(21) 



(22) 



Again the divergence of the integral is cut off by the distribution function. Now the peak is at non vanishing frequency 
and it disperses a to = £k 



A s ing (k, CJ) 



1 



"/go 
VAqo V |£ - w| 



(23) 



II. ROLE OF LONGITUDINAL QUANTUM PHASE FLUCTUATIONS 



In order to estimate the effect of longitudinal quantum phase fluctuations on the low energy fermion properties, we 
consider the following model of two dimensional nodal Dirac fermions coupled to gaussian phase fluctuations: 



S = S F + S e + XSi 



(24) 



Si 



(25) 
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Sg = dr 



d 2 x 



(d T 6) 2 + Ao (W) 2 



(26) 



Sl = I dT I !/Jj ^' ^r-v««i'.«i'. 



(27) 



Here tyj = (V'j.+.T: i s the Nambu spinor related to the jth pair of nodes at ±kj (hence the index ± of ip). A is 

a small parameter. The phase action Sg is taken to be strictly in 2d, rather than a layered system. This simplification 
does not change the results below, since the essential property of Sg is the linear 2d-like dispersion near the origin. 
As we discussed in the paper, the full layered dispersion is nearly linear near the origin (except for a small gap at 
k = due to the c-axis Josephson coupling, which is negligible at the temperatures of interest here). 

Si is the minimal coupling action of the fermions to the phase field, of which we keep only the most relevant term 
consisting of a current-current coupling. 

We will treat XSi perturbatively. The "engineering" scaling dimensions of the fields can be read off from the 
fixed point action with A = 0: [*f?j] = — 1, [9] = — |. Therefore the scaling dimension of Si is found to be [Si] = 
— i + O (A), i.e. it is irrelevant in the weak coupling limit. Physically, this means that the coupling to longitudinal 
phase fluctuations does not change the low energy spectrum of the fermions. In particular, we can estimate the 
quasiparticlc lifetime due to Si. To leading order, 1/t(E) ~ A 2 . In order to determine the energy dependence, we 
can perform an RG transformation that takes E to some fixed energy scale Eq > E. Under this transformation, 



A^ A' 



E 



-|+0(A) 



therefore 



1 



f „ n -1 + 0(A) 



Now, since 1/r has units of energy, we can scale back to get 



E 1 



r (E) E t (E ) 



E 



2-0(A) 



A 



(28) 



(29) 



(30) 



We conclude that for weak coupling 1/t(E) ~ E 2 , and the low energy quasiparticles are well defined. A direct 
evaluation of the leading diagram for 1/r (E) confirms this result. 

III. CALCULATION OF THE TYPICAL CURRENT FLUCTUATION 

The action for the transverse (vortex) part of the phase field in the 2d x — y model is the form of a Coulomb gas 
model, which can be mapped onto the Sine-Gordon model 



Ssg = / d 2 x 



^4^ (W)2 ~| COS[0(X)+ ^ (X)] 



(31) 



Here /x v is a vortex chemical potential term, that enables a calculation of vortex density correlations. In particular, 
the typical current fluctuation is given by 



^Hv.l(0))4E ( " (q) q f q)) 



(32) 



where n (q) is the Fourier transformed vortex density. The vortex density correlation function in real space can be 
expressed as: 



(q)M-q)> = ^/ |? e - q - x MxMo)) 

(n(x)n(O)) = -(2y) 2 (sin[0(x)]sin[0(O)])+2^(x)(cos[0(O)]) 



(33) 
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The last line is obtained by taking fi v derivatives of the Sine-Gordon partition function. Charge neutrality of the 
Coulomb gas model fixes / cPx (n (x) n (0)) = 0. Indeed, we see that otherwise (which is proportional to the total 
electostatic energy) diverges. 

In order to calculate Ij33|) in the gapped (disordered) phase, we use the self consistent harmonic approximation. 
The action (f3~Tj) is replaced by the quadratic action 



S SCHA = ^ / d 2 X 



(34) 

to is a variational parameter determined by minimizing the free energy of the system. The optimal value is to = 



j- (&yK) 2 v ^1 . The calculation of (|33|) with the action |34|) is strait forward, and yields 

(n(x)n(0)) = -(2 2/ ) 2 (^TO) 2 ^sinh[2 7 rA7(m|x|)]+^ c 2 ( 5(x) 
r fir, piqx 

/(m|x|) = / ? (35) 

^ ' U 2tt q 2 +TO 2 V ; 

Here A = 2y (cos [0 (xi)]). Instead of calculating A directly, we will adjust it so that (|3"5|) satisfies the exact charge 
neutrality sum rule / cPx (n (x) n (0)) = 0. In the high temperature (K — > 0) limit, p5|) reduces to the Debye-Hiickel 
form. 

Plugging ([33)) into (|3"2"|) and using J ^ ^JL — m ^ (where L is the system size, which is taken to infinity at the 
end of the calculation), we finally get 



Qo (T) = {Z c m?* K l°° drr In r sinh [2wKf (m&r)] (36) 

I 



